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ADAPTATION OF THE THEODORSEN THEORY TO THE REPRESENTATION 
OF AN AIRFOIL AS A COMBINATION OF A LIFTING LINE 
AND A THICKNESS DISTRIBUTION 

Raymond L. Barger 
Langley Research Center 

SUMMARY 

A representation of the Theodorsen airfoil theory as a combination of a lifting line 
and a thickness distribution is described. The approximations of thin -airfoil theory are 
avoided, since the full potential theory is used throughout. The theory provides a direct 
method for resolving an airfoil into a lifting line and a thickness distribution as well as a 
means of synthesizing thickness and lift components into a resultant airfoil and computing 
its aerodynamic characteristics. Specific applications of the technique are discussed. 

INTRODUCTION 

The thin-airfoil technique of representing an airfoil as a combination of a thickness 
distribution and a camber line is widely used because of the availability of a body of appro- 
priate data (ref. 1) and because the method facilitates design, inasmuch as it enables the 
designer to treat the lifting and thickness components separately. On the other hand, the 
thin-airfoil theory itself is vulnerable to criticism on several grounds regarding both the 
consistency of the theory and the accuracy of the procedures. For example, the procedure 
of superimposing a thickness distribution normal to the camber line involves a nonuniform 
stretching of the basic thickness profile. Furthermore, the simple addition of velocities 
due to thickness, camber, and angle of attack is not always reliable. An even more seri- 
ous difficulty is the failure of thin-airfoil theory to provide reasonable velocity values in 
the nose region. To improve the accuracy of the calculations in this region requires in- 
volved procedures. Finally, although the theory gives explicit equations for computing 
airfoil coordinates from thickness and camber distributions, it provides no analytic means 
for accomplishing the inverse task — that of resolving a given airfoil into its thickness 
and camber components. 

The present report describes a method for representing the Theodorsen airfoil 
theory as a combination of a lifting line and a thickness distribution. In this approach the 
approximations of thin-airfoil theory are avoided, since the full potential theory is used 
throughout. Within the framework of the full potential theory, the new method not only 



provides formulas for synthesizing an airfoil from a lifting line and a thickness distribu- 
tion, but it also provides analytic means for resolving a given airfoil into its thickness 
and lifting-line components. This latter capability is especially useful in certain airfoil 
design procedures. Both the synthesis and the analysis procedures are simple in concept 
and convenient in application. 


SYMBOLS 


An3n 


a 

Cp 

c 

R 


r 


V 


V 


x,y 


a 

“I 


Fourier coefficients 



pressure coefficient 
airfoil chord length 

radius of circle into which an airfoil is mapped by the Theodorsen 
transformation 

radial coordinate in near -circle plane 
undisturbed free -stream velocity 
local velocity 
?irfoil coordinates 
angle of attack 

+ ^te 

ideal angle of attack, - — 


Q!q angle of attack at zero lift 

e function relating angular coordinates of near-circle and exact-circle airfoil 

transformations 

value of e at airfoil nose, e{6) at 0 = 0 
e|.g value of e at airfoil trailing edge, e(0) at 9 = n 
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0 angular coordinate in near-circle plane 

(/) angular coordinate in exact -circle plane 

(/)* dummy variable of integration 

xj/ function relating radial coordinates of near-circle and exact-circle airfoil 

transformations 

average value of 

Subscripts: 

a antisymmetric 

s symmetric 

Primes indicate derivatives with respect to 6. 

BASIC CONCEPTS 
Potential Theory 

Before reviewing some thin-airfoil techniques, it is convenient first to state the 
basic formulas of the Theodorsen transformation theory for airfoil analysis. These for- 
mulas are useful in understanding the following discussion of thin-airfoil techniques, and 
they are required in the subsequent analysis section. 

The Theodorsen airfoil theory (ref. 2) involves the Joukowski transformation of the 
airfoil into a shape approximating a circle. In the transformed plane, this near circle is 
described by polar coordinates (r,0). The approximate circle is then transformed into 
an exact circle having coordinates The function i// and the constant xjy^ are 

defined by r = ae^ and R = ae^o. Theodorsen shows that is the average value 

of x{/, that is, 
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and that the functions 'P - 3^nd e = (p - 6 

1 p277 (b* - d) 

if/ - If. = ^ e cot ** ' '*' d.f>* 
o 2 tt Jq 2 


are related by the conjugate equations 

( 2 ) 

(3) 


The function \p is directly related to the airfoil coordinates in the physical plane by 


X = 2a cosh rp cos 6 
y = 2a sinh ip sin 9 


(4) 


In practice, xp is obtained directly as a function of 9 from equations (4) , and then the 
function e(9) is obtained by taking the conjugate of xj/. To be more precise, e{4>) 
should be obtained by iteration by means of the relation p = e + 9, and then e{9) should 
be calculated from e{9) = e|0(0)]. However, for practical airfoil shapes the additional 
precision gained by this step is slight (ref. 3), and so it is often omitted. The velocity 
distribution on the airfoil is computed from the equation 


V 


j"sin(Q! + e + 9) + sin^O!+ (1 + e’je'^o 


^(sinh^i// + sin20)(l + x{/'^) 


(5) 


Comments on Thin -Airfoil Methods 

The thin -airfoil theory, which affords a simple approximate procedure of synthesiz- 
ing and analyzing an airfoil by combining lifting-line and thickness elements, does not pro- 
vide a means for separating a given airfoil into its camber line and thickness distribution. 
Furthermore, the thin -airfoil calculation of the velocity, by a distribution of sources and 
sinks along the axis, gives a poor approximation, especially in the nose region. Attempts 
to improve the accuracy of the basic thin -airfoil theory have led to complicated proce- 
dures. (See sec. 4.5 of ref. 1, and ref. 4.) 

Another difficulty is that the thin-airfoil method of obtaining a profile by superim- 
posing a thickness distribution (fig. 1(a)) normal to a camber line, as illustrated in fig- 
ure 1(b), is an approximate procedure. The arc length of the camber line is longer than 
the chord of the original thickness profile, and therefore the thickness form undergoes a 
distortion, the stretching being greater for the upper surface than for the lower surface. 
The simpler thin-airfoil procedure, which attributes the lift to a mean line with the thick - 
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ness distribution superimposed on it in a vertical direction (fig. 1(c)), is an even poorer 
approximation. With this procedure, not only is the thickness distribution distorted, but 
the effective camber line is changed also. 

A theory that avoids some of the difficulties of the thin-airfoil methods is described 
in the following section. 


ANALYSIS AND APPLICATIONS 

Airfoil Synthesis From Thickness and Lift Components 

This section describes the theory and procedure for synthesizing an airfoil by com- 
bining an arbitrary thickness distribution with an arbitrary lifting line within the context 
of the full potential theory. The aerodynamic relationships between the resulting airfoil 
and its thickness and lifting-line components are also discussed. 

In equations (a) and (b) of reference 2, it is shown that (if 0 is substituted for (p, 
as previously discussed) e and xp can be expressed by 


oo 

-e = y -^fB ^ cos n6 - Ajj sin n0) 
rh \ ^ 

n=0 


OO 

COS n0 + sin nO) 

n^O 


From its definition, e must be a continuous function of 6, and therefore, 


e(0) = e(2n) s 


It follows from equation (6a) that 



de = 0 


(6a) 


(6b) 


(7) 


( 8 ) 


Now consider first the case for which e is antisymmetric about it (fig. 2(a)); 
that is, e(iT-0) = -e( 7 T+e). Then = 0 for all n, and therefore the conjugate function 
^ is symmetric about it (fig. 2(a)). Since x{/^ is constant, xpijs-d) = xp{Tr+6). 

Consequently cosh = cosh i//(ir+0) and sinh 4/{ir-9) = sinh i//(tt+ 0). Thus in equa- 

tions (4) for the airfoil coordinates, x(n+6) = x(n-0); that is, x is symmetric with 


5 



respect to ir. But because the sin 6 factor is an odd function, y is antisymmetric; 
that is, y(7T+0) = -yiir-O). In other words, since 9 = it corresponds to the trailing edge, 
the airfoil is symmetric, and so it represents a thickness distribution (fig. 2(b)) with the 
same pressure distribution for upper and lower surfaces (fig. 2(c)). The maximum thick- 
ness is controlled by the average value of ip. 

It is easily verified that for symmetric airfoils (thickness distributions) the condi- 
tions expressed by equations (7) and (8) are assured by the antisymmetry of the e-function 
together with the continuity conditions. Furthermore it is seen from figure 2(a) that, in 
this case, both e^^ and e|,g vanish. Thus the angle of attack at zero lift = -e^^ 

+ ^te 

and the ideal angle of attack cKj = are both zero, as expected for a symmetric 

airfoil. These are important aerodynamic parameters because e^-g is proportional to 
the lift at o: = 0 and because a-^ is useful in locating the bucket of the drag curve (see 
pt. Ill of ref. 5). When the converse situation, for which e is symmetric about n 
(fig. 3(a)), is considered, xp is antisymmetric about tt (fig. 3(a)), provided that xp^ = 0. 
In that case, x{ir+6) = x{-n-6) because cosh is an even function of its argument. Also 
y(7T+6) = y(7T-0) because the product of the two odd functions sinh p and sin 9 is even. 
In other words the upper and lower surfaces are identical, and the airfoil consists of a 
single line (fig. 3(b)). Its pressure distribution is shown in figure 3(c). 

For such a lifting line the e-function must satisfy 
C'Tr 

6(e) de = o 

-'o 


in order to satisfy equation (8). The magnitude of the lift at a = 0 is determined by 
®o “ "^te- Usually ej^ has a somewhat larger negative value than cKq, so that 
^te . 

Q!j = IS positive. 

Now consider the e-function that results from adding the e-function for the lifting 
line to that of the symmetric airfoil. If neither e-function is identically zero, the sum is 
neither symmetric nor antisymmetric about it. (See fig. 4(a).) Since both component 
e-functions are continuous and satisfy the conditions expressed by equations (7) and (8), 
the resultant e-function also satisfies these conditions, and so e corresponds to a real 
airfoil. Also, since e^g and ej.^ are both zero for the symmetric airfoil, these param- 
eters have the same values in the resulting e-function as for the lifting line. This means 
that the airfoil corresponding to the resultant e-function has the same angle of attack at 
zero lift and the same design lift coefficient as that of the lifting line. 
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In order to construct an airfoil corresponding to the e-function obtained by adding 
those of a lifting line and of a thickness distribution, it is necessary to have the corre- 
sponding ;// - function. This function can be obtained as the conjugate of e, or, if 
the ;//-f unctions for both component airfoils are known, it is obtained simply by adding the 
\p - \j/^ functions for the components. The parameter should be assigned the same 

value as that for the symmetric airfoil since the value of xj/^ for the lifting line is zero. 
Then the airfoil coordinates are determined by equations (4), and the velocity distribution 
is computed from equation (5). Results for the sample case are shown in figures 4(b) 
and 4(c). 

The constant is basically a thickness parameter, but this does not necessarily 

mean that the maximum thickness of the resultant airfoil will be exactly the same as that 
of the original thickness distribution, although it should be very nearly the same. How- 
ever, the important factor, e^o, in the velocity equation (5) will be exactly the same as 
that for the thickness distribution. This factor is not a function of position on the airfoil, 
and so it determines the level of the velocity curves independently of the details of the 
shape of the airfoil. 

Thus it is seen that this method of combining a lifting line and a thiclmess distribu- 
tion to synthesize an airfoil has a more precise aerodynamic significance than the thin- 
airfoil methods. The resultant airfoil has the same angle of attack at zero lift and ideal 
angle of attack as the lifting line and the same velocity -amplitude factor e^® as the 
original thickness profile. Its pitching moment is not exactly the same as that of the lift- 
ing line, because the pitching moment is not entirely independent of the thickness distribu- 
tion according to the full potential theory. The procedure may be summarized as follows: 

(1) Determine the ^//-functions for the thickness and lifting-line components from 

their coordinates. Compute the average value of xp for the thickness distribution 

and then compute the conjugate e-functions from equation (2). 

(2) Add the e-functions and the ^//-functions to obtain those for the resultant airfoil. 

(3) Obtain the resultant airfoil coordinates from equations (4) and the velocity dis- 
tribution for this airfoil from equation (5). 

If one intends to maintain a catalog of thickness distributions and lifting lines, he 
should also maintain a file of the corresponding e- and -functions. The simple addition 
of a pair of these functions yields the resultant transformation functions from which both 
the coordinates and the velocity distribution for this airfoil can be obtained. 

Airfoil Resolution Into Thickness and Lift Components 

In order to find the lifting line and thickness distribution corresponding to a specified 
arbitrary airfoil (fig. 5), it is necessary to determine its transformation functions and to 
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resolve them into their symmetric and antisymmetric parts (fig, 6). The formulas (for 
the e-function) are 

e(0) + e(27T-e) 

= 2 

_ e(g) - e (277-0) 
a 2 

It is readily verified that es(7T-0) = es(7r+0) and e„{T7-9) = -e„{TT+0). 

H. a. 

The first of these functions, eoj corresponds to the lifting line, and e corre- 
spends to the thickness distribution. Conversely, when \p - ip ^ is resolved into its 
components, the symmetric part corresponds to the thickness distribution, and is 
set equal to that of the original airfoil. For the antisymmetric part of xp - ip^, corre- 
sponding to the lifting line, xp^ = 0. The xp- and e-functions for each of these components 
can then be used to compute the profiles from equations (4) (figs. 5(b) and 5(c)) and the 
velocity distributions from equation (5) (fig. 7). 

Application of the theory to altering certain properties of an airfoil can be demon- 
strated by a rather extreme example. The airfoil of figure 5(a) is taken as the basic 
shape. Its lift at zero angle of attack is increased by 50 percent by adding to the e-function 
of its lifting line (fig. 6(c)) a bilinear function that increases from -0.5e|.g at 0 = 0 to 
+0.5e^g at 0 = 77 and then decreases linearly to -0.5ej.g at 0 = 2t7. The thickness of the 
airfoil is increased from 10.7 percent to 15 percent by changing the value of xp^ for the 
thickness distribution (fig. 6(b)) from 0.10 to 0.137. The new e- and i//-functions are then 
combined, and the resulting airfoil and its pressure distribution are shown in figure 8. 

If it is desired to change the lifting-line component of an airfoil in order to redis- 
tribute the loading, the calculations should be carried out at the ideal angle of attack, at 
which the leading-edge velocity is not infinite for the lifting line. For example, the pres- 
sure distributions shown in figure 9(a) (for the airfoil) and in figure 9(b) (for its lifting 
line), which are computed at aj, may be compared with those at a = 0 for the same 
airfoil and for its lifting line (figs. 4(c) and 3(c), respectively). 

CONCLUDING REMARKS 

A representation of the Theodorsen airfoil theory as a combination of a lifting line 
and a thickness distribution is described. The approximations of thin-airfoil theory are 
avoided, since the full potential theory is used throughout. The theory provides a direct 
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method for resolving an airfoil into a lifting line and a thickness distribution as well 
as a means of synthesizing thickness and lift components into a resultant airfoil and 
computing its aerodynamic characteristics. Some specific applications of the technique 
are discussed. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, Va. 23665 
December 2, 1975 
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(a) Thickness distribution. 



(b) Lifting line used as camber line with thickness 
superimposed normal to the line. 



(c) Lifting line used as mean line with thickness 
superimposed vertically. 

Figure 1.- Thin-airfoil procedures for combining a thickness 
distribution with a lifting line. 
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(a) Transformation functions. 




x/c 

(c) Pressure distribution, a = 0, 

Figure 2.~ Calculations for symmetric airfoil (thickness distribution). 





(a) Transformation functions. 


(b) Lifting -line airfoil shape. 



(c) Pressure distribution, a = 0. 
Figure 3.- Calculations for lifting -line airfoil. 
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(a) Transformation functions obtained by adding functions of figure 2(a) 
to corresponding functions of figure 3(a). 



(b) Airfoil profile. 



(c) Pressure distribution, a = 0. 

Figure 4.- Theoretical results for airfoil synthesized from thickness 
distribution of figure 2 and lifting line of figure 3. 
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(a) Given airfoil. 



(b) Thickness distribution. 



(c) Lifting line. 

Figure 5.- Resolution of an airfoil into thickness and 
lifting-line components. 
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(a) Functions corresponding to given airfoil. 



0 


(b) Symmetric part of \{y and antisymmetric part 
of e (thickness component). 



0 


(c) Symmetric part of e and antisymmetric part 
of (lifting “line component). 

Figure 6.- Transformation functions corresponding respectively to 
airfoils of figures 5(a), 5(b), and 5(c). 




(a) For given airfoil. 




(c) For lifting line. 


Figure 7. - Pressure distributions corresponding respectively 
to airfoil shapes of figures 5(a), 5(b), and 5(c). 







NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
WASHINGTON. D.C. 20546 


OFFICIAL BUSINESS 

PENALTY FOR PRIVATE USE $300 SPECIAL FOURTH-CLASS RATE 

BOOK 

V 

t /•' 


POSTAGE AND FEES PAID 
NATIONAL AERONAUTICS AND 
SPACE ADMINISTRATION 
451 



rjg 001 Cl u a 751212 soososds 

DEPT OF THE AIR FORCE 
AF WEAPONS LABORATORY 
ATTN: TECHNICAL LIBRARY (SUL) 
KIRTLAND AFB NM 87117 


POSTMASTER ; 


If Undeliverable (Section 158 
Postal Manual) Do Not Return 


^'T/?e aeronautical and space activities of the United States shall he 
conducted so as to contribute to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof!* 

— National Aeronautics and Space Act of 1958 


NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 


TECHNICAL REPORTS: Scientific and 

technical information considered important, 
complete, and a lasting contribution to existing 
knowledge. 

TECHNICAL NOTES: Information less broad 
in scope but nevertheless of importance as a 
contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: 

Information receiving limited distribution 
because of preliminary data, security classifica- 
tion, or other reasons. Also includes conference 
proceedings with either limited ^ unlimited 
distribution. 

CONTRACTOR REPORTS: Scientific and 
technical information generated under a NASA 
contract or grant and considered an important 
contribution to existing knowledge. 


TECHNICAL TRANSLATIONS: Information 
published in a foreign language considered 
to merit NASA distribution in English. 

SPECIAL PUBLICATIONS: Information 
derived from or of value to NASA activities. 
Publications include final reports of major 
projects, monographs, data compilations, 
handbooks, sourcebooks, and special 
bibliographies. 

TECHNOLOGY UTILIZATION 
PUBLICATIONS: Information on technology 
used by NASA that may be of particular 
interest in commercial and other non-aerospace 
applications. Publications include Tech Briefs, 
Technology Utilization Reports and 
Technology Surveys. 


Details on the availability of these publications may be obtained from: 

SCIENTIFIC AND TECHNICAL INFORMATION OFFICE 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington, D.C. 20546 



